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We present explicit examples to show that the 'compatibility criterion' [recently ob- 
tained by us towards providing equilibrium configurations compatible with the struc- 
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ture of general relativity] which states that: for a given value of o~[= (Pq/Eq) = the 
ratio of central pressure to central energy-density], the compactness ratio u[= (M/R), 
' where M is the total mass and R is the radius of the configuration] of any static config- 
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uration cannot exceed the compactness ratio, u^, of the homogeneous density sphere 
(that is, u < Uh), is capable of providing a necessary and sufficient condition for 
any regular configuration to be compatible with the state of hydrostatic equilibrium. 
This conclusion is drawn on the basis of the finding that the M — R relation gives the 



necessary and sufficient condition for dynamical stability of equilibrium configura- 



tions only when the compatibility criterion for these configurations is appropriately 
satisfied. In this regard, we construct an appropriate sequence composed of core- 
envelope models on the basis of compatibility criterion, such that each member of 
this sequence satisfies the extreme case of causality condition v = c = 1 at the centre. 
The maximum stable value of u ~ 0.3389 (which occurs for the model corresponding 
to the maximum value of mass in the mass-radius relation) and the corresponding 
central value of the local adiabatic index, (ri)o ~ 2.5911, of this model are found 
fully consistent with those of the corresponding absolute values, it m ax 

< 0.3406, and 

(ri)o < 2.5946, which impose strong constraints on these parameters of such models. 
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In addition to this example, we also study dynamical stability of pure adiabatic poly- 
tropic configurations on the basis of variational method for the choice of the 'trial 
function', £ = re u / 4 , as well as the mass-central density relation, since the compati- 
bility criterion is appropriately satisfied for these models. The results of this example 
provide additional proof in favour of the statement regarding compatibility criterion 
mentioned above. Together with other results, this study also confirms the previous 
claim that just the choice of the 'trial function', £ = re u ^ 4 , is capable of providing 
the necessary and sufficient condition for dynamical stability of a mass on the basis 
of variational method. Obviously, the upper bound on compactness ratio of neutron 
stars, u = 0.3389, which belongs to two-density model studied here, turns out to 
be much stronger than the corresponding 'absolute' upper bound mentioned in the 
literature. 

PACS Nos.: 04.20.Jd; 0440.Dg; 97. 60. J d. 
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1. INTRODUCTION 



Einstein's field equations for static and spherically symmetric mass distribution were first 
solved by Schwarzschild (1916). The first solution describes the geometry of the space- 
time exterior to a prefect fluid sphere in hydrostatic equilibrium. While the other, known 
as interior Schwarzschild solution, corresponds to the interior geometry of a fluid sphere 
of constant (homogeneous) energy-density, E. The importance of these two solutions in 
general relativity is well known. The interior Schwarzschild solution provides two very 
important features towards obtaining configurations in hydrostatic equilibrium, compatible 
with general relativity, namely - (i) It gives an absolute upper limit on compactness ratio, 
u{= M/R, mass to size ratio of the entire configuration in geometrized units) < (4/9) for 
any static and spherical configuration (belonging to arbitrary density profiles, provided the 
density does not increase outwards) in hydrostatic equilibrium (Buchdahl 1959; Weinberg 
1972), and (ii) For an assigned value of the compactness ratio, u, and radius R (or mass M), 
the minimum central pressure, Po, corresponds to the homogeneous density solution (see, 
e.g., Weinberg 1972). 

Despite the non linear and coupled differential equations, various exact solutions of the 
field equations for static and spherically symmetric metric are available in the literature 
(see, e. g., Kramer et al 1980) which may be used to obtain various physical properties 
of spherical and static compact object (provided they are physically realistic). Knutsen 
(1988; 1989) examined physical properties of the various exact solutions and found that 
these solutions correspond to nice physical properties and also remain stable against small 
radial pulsations upto certain values of u. 

Another way to explore the physical properties of compact objects like neutron star, 
one may expect to have some physically viable equation of state (EOS). However, for such 
objects the equations of state (EOSs) are not well known [empirically] because of the lack 
of knowledge of nuclear interactions beyond the density ~ 10 14 gcm -3 (Dolan 1992), and 
the only way to obtain EOSs far beyond this density range is extrapolation. Various such 
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extrapolated equations are available in the literature (Arnett & Bowers 1977). As a way 
out, one can impose some restrictions upon the known physical quantities, such that, the 
speed of sound inside the configuration, v[= \J\dP '/ 'dE)\, does not exceed the speed of light 
in vacuum, i.e., v < c = 1 (in geometrized units), and obtain various physical properties, like 
upper bound on stable neutron star masses (Rhoades & Ruffini 1974; Brecher & Caporaso 
1976; Hartle 1978; Friedman & Ipser 1987). Haensel and Zdunik (1989) have shown that 
the only EOS which can describe a submillisecond pulsar and the static mass of 1.442M 
simultaneously, corresponds to the EOS, (dP/dE) = 1, however, they emphasized that this 
EOS represents an 'abnormal' state of matter in the sense that pressure vanishes at densities 
of the order of nuclear density or even higher. 

We have recently proposed a core-envelope model with stiffest EOS, (dP/dE) = 1, 
forming the core and a polytropic equation with constant adiabatic index r\ = {dlnP/ dlnp) 
[where P is the pressure and p represents the rest-mass density] describing the envelope, such 
that P, E, both of the metric parameters (z/, and A), their first derivatives, and the speed of 
sound are continuous at the core-envelope boundary and at the exterior boundary (surface) 
of the structure (Negi & Durgapal 2000). The other remarkable feature of this core-envelope 
model is that not only the 'abnormalities' (in the sense discussed in the literature, see, e.g. 
Lee 1975; Haensel & Zdunik 1989) disappear, the maximum value of u = 0.3574 for the stable 
configuration turns out to be as large as that obtained by using the EOS, (dP/dE) = 1, 
alone (see, e. g., Haensel & Zdunik 1989). The maximum stable mass of neutron star based 
upon this model (by using the maximum value of u = 0.3574 for stable configuration) turns 
out to be 7.944M , if the (average) density of the configuration is constrained by fastest 
rotating pulsar, with rotation period, P rot = 1.558 ms, known to date. The model gives 
pulsationally stable configurations with compactness ratio u > (1/3), which are important 
to study Ultra-Compact Objects (UCOs) [see, e. g., Negi & Durgapal 1999a, b; 2000; and 
references therein]. 

Recently, by using property (ii) of the homogeneous density sphere as mentioned above, 
we have connected the compactness ratio, u, of any static and spherical configuration with 
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the corresponding ratio of central pressure to central energy-density a[= (P /E )] and 
worked out an important criterion which concludes that for a given value of a, the maximum 
value of compactness ratio, u(= Uh), should always correspond to the homogeneous density 
sphere (Negi & Durgapal 2001). 

An examination of this criterion on some well known exact solutions and EOSs indicated 
that this criterion, in fact, is fulfilled only by two types of configurations corresponding to 
a single EOS or density variation: (i) the regular (positive finite density at the origin which 
decreases monotonically outwards) configurations which correspond to a vanishing density 
at the surface together with pressure [so called, the gravitationally-bound structures] (Negi 
& Durgapal 2001; Negi 2004b; Negi 2006), and (ii) the structures which correspond to a 
non-vanishing surface density but exhibit singularities at the centre, in the sense that both 
pressure and density become infinity as r — > [so called, the self-bound singular structures] 
(Negi 2004b; Negi 2006). On the other hand, it is seen that the EOSs or analytic solutions, 
corresponding to a non-zero finite, surface density (that is, the pressure vanishes at finite 
surface density, and so called the self-bound regular structures), in fact, do not fulfill this 
criterion (Negi & Durgapal 2001; Negi 2004b; Negi 2006). We have shown this inconsistency 
particularly for the EOS, (dP/dE) = 1 (as it represents the most successful EOS to obtain 
the various extreme characteristics of neutron stars as discussed above). 

In addition to the self-bound regular configurations corresponding to a single density 
variation, the compatibility criterion may not be satisfied by various two-density, or multiple- 
density, regular, gravitationally-bound structures. Such structures are widely discussed in 
the literature, particularly, for determining the upper bound on neutron star (NS) masses. 
In this connection, we would consider the core-envelope model proposed by Negi & Durgapal 
(2000). 

The reason(s) behind non-fulfillment of the criterion obtained in the study of Negi and 
Durgapal (2001) by various exact self-bound regular solutions and EOSs, as well as the 
two-density, gravitationally-bound, regular structures, and their further implications are 
discussed in the following sections. 
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2. Compatibility criterion: the necessary and sufficient condition for hydrostatic equi- 
librium of any static spherical configuration 



In order to bring things together regarding the compatibility of regular structures mentioned 
above, we follow Negi and Durgapal (2001) by assuming a homogeneous sphere of uniform 
energy-density, E. The equations for isotropic pressure P, and uniform energy-density E, 
can be written in terms of compactness ratio, u, and the radial coordinate measured in units 
of configuration size, y{= r/R) as 



Let us consider a regular variable density sphere (with some given EOS or analytic 
solution) with central energy-density Po and central pressure Po , corresponding to the 
compactness ratio u = u v . 

Now, we can always construct a homogeneous density sphere with the same value of 
the compactness ratio u v , and energy-density E , because if Pot corresponds to the central 
pressure of this sphere, the ratio <J^(= Poh/E ) depends only upon the assigned value of the 
compactness ratio u v . And, P h is given by 



8nER 2 = 6u. 



(1) 



8irPR 2 = 6u 



(1 -2m/ 2 ) 1 / 2 - (1 -2m) 1 / 2 
3(1 - 2u) 1/2 - (1 - 2uy 2 ) 1/2 



(2) 



Po. = (Qu/87tR 2 )\(1 - (1 - 2 M ) 1 / 2 )/(3(l - 2m) 1 / 2 - 1)" . 



(3) 



Now, according to property (ii) of homogeneous density sphere, we may write 



Po > P 



Oil 



(4) 



or, 



(Po/Po) > (Pofc/Po)- 



(5) 



Hence for a given value of u{= u v ), we obtain 



(6) 
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where a v is defined as the ratio, (Pq/Eq). 

Now, varying the compactness ratio, u v , for the homogeneous density sphere from u v to 
Uh (say), such that, we should have 

cr v = <?h- (7) 
For u = Uh , the value of ah would become 



((l-2« fc ) 1/2 -l)/(l-3(l-2u fc ) 1/2 )l- (8) 



Substituting Eq. (8) with the help of Eq. (7) into Eq. (6), we get 

(9) 



;i - 2« fc ) 1/2 - 1) i > r K(i - 2« w ) 1 /2 _ i> 



(1-3(1 -2« h )V2) J - L(i_ 3(1 -2« W )V2) 
Thus, it is clear from Eq. (9) that 

Uh > u v ( for an assigned value of a) . (10) 

That is, for an assigned value of the ratio of central pressure to central energy-density 
a(= a v ), the compactness ratio of homogeneous density distribution, u(= Uh) should always 
be larger than or equal to the compactness ratio u(= u v ) of any regular solution*, compatible 
with the structure of general relativity. Or, in other words, for an assigned value of the 
compactness ratio, u, the minimum value of the ratio of central pressure to central energy- 
density, a, corresponds to the homogeneous density sphere. 

In the light of Eq. (10), let us assign the same value M for the total mass corresponding 
to various regular configurations in hydrostatic equilibrium. If we denote the density of the 
homogeneous sphere by Eh , we can write 



*Notice that this rinding is also true for self-bound singular solutions because the ratio of (infinite) 
central pressure to density turns to be finite [examples of such solutions are well represented by 
Tolman's type V and VI solutions (Tolman 1939)]. Hence, the notion 'any regular' solution may 
be replaced by 'any static' solution [of course, with the requirement that the density decreases 
monotonically outwards from the centre] . 
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E h = 3M/(4 7 r J R fe 3 ) 



(11) 



where Rh denotes the radius of the homogeneous density sphere. If R v represents the radius 
of any other regular sphere for the same mass M, the average density E v of this configuration 
would correspond to 



That is, for an assign value of a the average energy- density of any regular configuration, E v , 
should always be less than or equal to the density, E h , of the homogeneous density sphere 
for the same mass M. 

We point out that the regular configurations corresponding to a single exact solution, 
or EOS with a finite central and non- vanishing surface density, in fact, do not fulfill the 
definition of this 'actual' total mass, M, which appears in the exterior Schwarzschild solution 
[this definition asserts that, being the coordinate mass, the particular 'type' of density 
variation considered for it should remain 'unknown' to an external observe and this is possible 
only when the mass depends either upon the central density, or upon the surface density, 
and in any case, not upon both of them (Negi 2004b; Negi 2006). It follows, therefore, that 
the central density should be independent of the surface density or vice-versa, according 
to the density distribution assigned for the mass]. And the so called self-bound regular 
structures, in fact, violate this requirement as they correspond to a surface density which 
always depends upon the central density and vice-versa. Thus, the main findings of the 
study regarding this criterion can be summarized in the following manner: 

(a) The gravitationally-bound regular configuration and self-bound singular structures, 
described by a single EOS or exact solution, fulfill the definition of the total mass, M, appears 
in the exterior Schwarzschild solution, hence the condition of hydrostatic equilibrium is 
naturally satisfied by these structures [this finding is fully consistent with the 'compatibility 



E v = 3M/(4tt j R„ 3 ). 



(12) 



Eq. (10) indicates that R v > R h . By the use of Eqs. (11) and (12) we find that 



E v < E h . 



(13) 
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criterion' , because for all (possible) values of a, the condition u < Uh is fully satisfied by 
these configurations]. 

(b) The self-bound regular configurations, described by a single EOS or exact solution, 
can not fulfill the definition of the total mass, M, appears in the exterior Schwarzschild 
solution, as a result, the state of hydrostatic equilibrium can not be satisfied by them [this 
finding is also fully consistent with the 'compatibility criterion' , because such configurations 
correspond to the condition, u > Uh (for all possible values of a)]. 

(c) The only regular configuration which can exist under the category (b) mentioned 
above is described by the homogeneous density distribution. 

Note that the two-density or multiple-density models (that is, the structures governed 
by two or more EOSs assigned for different regions with appropriate matching conditions 
at the core-envelope boundaries) of both of the categories, (a) and (b) described above 
(such that the definition of the mass M mentioned above is appropriately satisfied) are 
quite possible, however, as we will show in the present paper that the fulfillment of the 
definition of the mass 'M' for any two-density model represents only a necessary condition for 
hydrostatic equilibrium, because the 'compatibility criterion' may not be satisfied by them. 
As we have noted earlier that the necessary condition for hydrostatic equilibrium (that is, 
the fulfillment of the definition of the mass M) put forward by the exterior Schwarzschild 
solution is also sufficient for a single EOS or exact solution assigned for the mass, because 
this fact is also supported by the 'compatibility criterion'. It follows therefore that the 
'compatibility criterion' is capable of ensuring a sufficient and necessary condition for any 
structure (including two-density or multiple-density distribution) in the state of hydrostatic 
equilibrium. 

To elaborate this statement more clearly, let us consider the core-envelope model dis- 
cussed by Negi and Durgapal (2000). The core of this model is described by an EOS which 
belongs to the category '(b)' mentioned above, and the matching of various parameters 
at the core-envelope boundary is assured by characterizing an envelope which belongs to 
the category '(a)' EOS. That is, 'overall' the model describes a gravitationally bound two- 
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density structure [of category 'a' mentioned above], such that the necessary condition for 
hydrostatic equilibrium put forward by exterior Schwarzschild solution at the surface of the 
configuration (that is, the mass, M, depends only upon the central density, meaning thereby 
that the definition of mass is appropriately satisfied), even then, the compatibility criterion 
for hydrostatic equilibrium (Negi & Durgapal 2001) turns out to be unsatisfied for this 
model (as shown under section 4 of the present study). Thus, it follows that the fulfillment 
of necessary condition for hydrostatic equilibrium at the surface, and the achievement of 
proper matching conditions at the core-envelope boundary are not sufficient to assure the 
condition of hydrostatic equilibrium for any two-density structure. However, the fulfillment 
of compatibility criterion alone could provides a necessary and sufficient condition for any 
regular configuration (including two-density structures) to be consistent with the state of 
hydrostatic equilibrium. 

In order to verify this statement, we would re- investigate the core-envelope model put 
forward by Negi and Durgapal (2000), based upon the said compatibility criterion for hydro- 
static equilibrium, such that for each (possible) assigned value of a, the compactness ratio 
of the whole configuration, u, remains less than or equal to the compactness ratio, Uh, of the 
corresponding sphere of homogeneous density distribution. Such an investigation is possible, 
because we can re-adjust the boundary, r&, of the core-envelope model in such a manner that 
for an assigned value of a, the 'average density', £' a „(say), of the whole configuration always 
remains less than or equal to the density, E h , of the homogeneous density sphere for same 
mass M. Thus, this criterion should be fulfilled by any regular configuration specified by 
a single density distribution, a core-envelope model, a core-mantle-envelope model, or any 
other complicated distribution of matter composed of various regions inside the configura- 
tion, in order to fulfill the state of hydrostatic equilibrium. This statement is verified on the 
basis of dynamical stability of some regular configurations, consistent with the compatibility 
criterion, in the following sections, and the results which are summarized as Theorem 2 and 
its subsequent corollaries in the following section, may be stated in the general form as the 
following theorem. 
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Theorem 1: The necessary and sufficient condition for hydrostatic equilibrium of any 
static* and spherical configuration is that for an assigned value of the ratio of central pressure 
to central energy-density, the compactness ratio u(= M/R) of the said configuration should 
not exceed the compactness ratio Uh of the corresponding sphere of homogeneous density 
distribution. 

3. Necessary and sufficient condition for hydrostatic equilibrium and dynamical sta- 
bility of regular configuration 

The absolute values are obtained by using a 'compressible' sphere of homogeneous energy- 
density (Negi 2004a), such that the following relation holds good for a constant I\ 

r\p dp 

P + E ~ dE' 

And the adiabatic speed of sound, v = ^dP/dE), becomes finite inside this configura- 
tion for a finite (constant) IV In order to satisfy the extreme case of causality condition 
v — c — 1 at the centre of this sphere, we obtain (P /E ) = 0.6271, which correspond to 
a u value = 0.3406, and the (critical) constant I\ = (ri) = 2.5946 respectively, for the 
dynamically stable configuration. This value of u(= 0.3406) represents an absolute upper 
bound, consistent with causality and dynamical stability, since it follows from the compati- 
bility criterion that for this maximum value of u, the corresponding value of (P /E ) of any 
regular configuration can not be less than 0.6271. Now, this result may be generalized for 
the sequences, composed of NS models such that every member of this sequence satisfies 
(dP/dE)o = 1 (here and elsewhere in the paper, the subscript '0' represent the value of the 
corresponding quantity at the centre), in the following manner that the maximum stable 

*the notion 'any static' instead of 'any regular' is used here in the general sense, since the dy- 
namical stability of singular solutions (which may also satisfy the compatibility criterion) does not 
correspond to any solution. 
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value of u (corresponding to the case of first maxima among masses in the M — R relation) 
and the corresponding central value of local (ri)[= (r^o] of such sequences must satisfy the 
inequalities -Umax < 0.3406, and (ri) < 2.5946 respectively, in order to ensure the necessary 
and sufficient condition for dynamical stability of a mass. Since these absolute values are 
obtained by using the 'trial function', £ = re"/ 4 , which is able to provide the necessary 
and sufficient condition for dynamical stability * in the variational method (Chandrasekhar 
1964a, b). 

The equilibrium sequences of the type mentioned here, in fact, are widely discussed 
in the literature, but not on the basis of compatibility criterion. The core-envelope models 
presented by Negi & Durgapal (2000) also represent an equilibrium sequence of this type. We 
re-constructed this sequence on the basis of compatibility criterion for the first time (for other 
example of such a sequence, consistent with the compatibility criterion, see, e.g. Negi 2005), 
and it is seen that the maximum stable value of compactness ratio and the corresponding 
central value of local (r x )[= (r^o] of such sequences are found fully consistent with those 
of the values obtained mentioned above (section 4; see also Negi 2005). It follows, therefore, 
that the M — R relation (or the mass-central density relation) provides the necessary and 
sufficient condition for dynamical stability of a mass only when the compatibility criterion 

^although, the variational method gives only a sufficient condition for the dynamical stability of 
a mass, since the results depend somewhat on the choice of a particular 'trial function'. However, 
we have shown that the choice of the particular trial function, £ = re^ 4 , is capable of providing 
the most rigorous results among the various trial functions (Negi & Durgapal 1999b; Negi 2004a), 
and because of this reason it could provide the necessary and sufficient condition for dynamical 
stability. It would not be out of place here to point out that the variational method gives precise 
results only for the configurations which correspond to a smooth variation of density from centre 
to the surface (see, e.g. Bardeen et. al 1966), and for this reason we did not consider this method 
for analyzing the stability of core-envelope models considered here. 
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for the equilibrium configuration is also satisfied. Or, in other words, the compatibility 
criterion is able to provide the necessary and sufficient condition for hydrostatic equilibrium 
of regular configurations. 

In order to verify the last claim, irrespective of the particular type of core-envelope 
models considered in the present study, we would further consider (section 5) the dynamical 
stability of pure polytropic configurations (P = Kp Tl ) on the basis of variational method for 
the choice of trial function f = re v/ \ as well as the mass-central density relation for some 
assigned values of constant I\, since the polytropic configurations appropriately satisfy the 
'compatibility criterion' (Negi & Durgapal 2001). 

The results which would follow from the study of sections 4 and 5 respectively, may be 
summarized as theorem 2 and its subsequent corollaries in the following form 

Theorem 2: The mass-radius (or, mass-central density) relation for an equilibrium se- 
quence of regular configurations provides the necessary and sufficient condition for dynami- 
cal stability only when the equilibrium sequence it-self satisfies the necessary and sufficient 
condition of hydrostatic equilibrium (theorem 1). 

Corollary 1 to theorem 2: If an equilibrium sequence, composed of neutron star models in 
such a manner that every member of this sequence satisfies the extreme causality condition 
v — c — 1 at the centre, then the maximum value of compactness ratio (corresponding to 
the case of first maxima among masses in the M — R relation) and the corresponding central 
value of local adiabatic index (r^o are constrained by the inequalities, u max < 0.3406 and 
(Ti)o < 2.5946 respectively. 

Corollary 2 to theorem 2: For regular configurations, corresponding to enough smooth 
density variations such that the variational method could be used, the variational method 
could provide the necessary and sufficient condition of dynamical stability just for the choice 
of a particular trial function £ = re"'' 4 . 
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4. Hydrostatic equilibrium and dynamical stability of core-envelope models 



The metric for spherically symmetric and static configurations can be written in the following 
form 

ds 2 = e v dt 2 - e x dr 2 - r 2 d6 2 - r 2 sin 2 ^0 2 , (14) 

where v and A are functions of r alone. Recalling that we are using 'geometrized units', 
the Oppenheimer-Volkoff (O-V) equations (Oppenheimer & Volkoff 1939), resulting from 
Einstein's field equations, for systems with isotropic pressure P and energy-density E can 
be written as 

P' = -(P + E)[AnPr 3 + m]/r(r -2m) (15) 
v'/2 = -P'/(P + E) (16) 
m'{r) = A-kEt 2 ; (17) 

where m(r) is the mass contained within the radius r, and the prime denotes radial derivative. 

The core-envelope model (Negi & Durgapal 2000) consists of a core with most stiff EOS 
in the region < r < b, and an envelope with a polytropic EOS in the region b < r < R, 
given as 

(i) The core: < r < b 

For the models of neutron stars considered here, we have chosen the core of most stiff 
material as 

P = (E- E a ) (18) 

where E s is the value of density at the surface of the configuration, where pressure vanishes. 

(ii) The envelope: b < r < R 

The envelope of this model is given by the equation of state 

P = Kp Fl (19) 
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or 

(E-p) = P/(T l -l). 

where if is a constant to be worked out by the matching of various variables at the 
core-envelope boundary and p and r x represent respectively, the rest-mass density and the 
(constant) adiabatic index as defined earlier. 

At the boundary, r = b, the continuity of P(— P b ), E(= E b ), and r(= r b ) require 

^ = p 6 /(E b -[p 5 /(r 1 -i)]) ri (20) 

where Ti is given by (see, e.g., Tooper 1965) 
T 1 = [(P + E)/P](dP/dE). 

The continuity of (dP/dE), at the boundary gives 

r\ = 1 + (E b /p b ). (21) 

Thus, the continuity of P,E,v,\, and (dP/dE) at the core-envelope boundary is ensured, 
for the static and spherically symmetric configuration. 

The coupled Eqs. (15), (16), (17), are solved along with Eqs.(18) and (19) for the 
boundary conditions (20) and (21) [at the core-envelope boundary, r = b], and the boundary 
conditions, P = E = , m(r = R) = M, e v = e~ x = (1 - 2M/R) = (1 - 2u) at the external 
boundary, r = R. 

For the sake of numerical simplification, we assign the central density, E — 1. It is seen 
that the degree of softness of the envelope is restricted by the inequality, (P b /E b ) > 0.014. 
For the minimum value of (P b /E b ) = 0.014, we obtain various quantities, such as, core 
mass, M b , core radius, r b , density at the core-envelope boundary, E b , total mass, M, and 
the corresponding radius, R, of the configuration in dimensionless form. Some of these 
quantities are shown in Table 1 for various assigned values of the central pressure to density 
ratio, (Pq/Eq). 

To determine the stability of the models given in Table 1, we need to draw the mass- 
radius diagram for the structures. For this purpose, we have normalized the boundary 
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density, E b = 2 x 10 14 g cm" 3 , and obtained the mass-radius diagram as shown in Fig. 1 of 
Negi and Durgapal (2000) [Notice that the value of E b chosen in this way (and hence also the 
mass and the radius obtained in conventional units as shown in Fig. 1) is purely arbitrary. 
These values have nothing to do with the actual maximum mass and the corresponding 
radius of the stable neutron star obtained in the present paper.]. The maximum stable 
value of u of the whole configuration is obtained as 0.3574. For this maximum value of u, 
the binding energy per baryon, a r [= (M r — M)/M r , where M r is the rest-mass (Zeldovich 
& Novikov 1978) of the configuration] also approaches maximum (= 0.2441) as shown in 
Table 1. Although, the corresponding P /-^o(— 0.704) value (or (ri) ^ 2.4204 value) is 
consistent with the corresponding absolute value (P /E is larger than 0.6271 or (r^o is less 
than 2.5946), the configuration is not consistent with the corollary 1 of theorem 2, since the 
maximum value of u ~ 0.3574 is inconsistent with the absolute upper bound onu = 0.3406. 
It follows, therefore, that the M — R relation does not provide the necessary and sufficient 
condition for dynamical stability of equilibrium masses, since these models are not consistent 
with the compatibility criterion, which is evident from Table 1. As the first column of Table 
1 corresponds to compactness ratio, Uh, of homogeneous density sphere as calculated from 
Eqs. (1) and (2) for various assigned values of a shown in Table 1. Column seventh of this 
table represents the compactness ratio u of the whole configuration for the same values of 
a. Comparing column one and seventh, we find that for each assigned value of a, u > Uh, 
meaning thereby that the model is inconsistent with the compatibility criterion. 

To make the model consistent with the structure of general relativity, we re-investigate 
this model based upon the compatibility criterion (section 2) by solving the coupled Eqs. 
(15), (16), and (17) together with Eqs. (18) and (19) for the boundary conditions (20) 
and (21) respectively. For numerical simplicity, we assign the central energy- density of the 
configuration, E — 1. The order of numerical precision is set precisely following the specific 
nature of EOSs for the core and envelope regions respectively. The ratio, (Pb/E b ), at the 
core-envelope boundary is so adjusted that for each and every (possible) assigned value of 
<7, the compactness ratio, u, of the whole configuration always turns out to be less than or 
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equal to the compactness ratio, Uh, of the homogeneous density sphere for same values of 
a. The results obtained in this regard are shown in Table 2. It is seen that to meet the 
requirement set up by compatibility criterion, the minimum value of (Pb/Eb) reaches about 
2.9201 x 1CT 1 . 

To investigate the stability of the models which are now compatible with the structure of 
general relativity and causality (Table 2), we draw the mass-radius diagram for the models 
by normalizing the boundary density, Eb = 2 x 10 14 gcm~ 3 , as shown in Fig. 1 [notice 
that the use of normalizing density, Eb, as mentioned in the previous case also, is purely 
arbitrary and its purpose is only to determine the maximum value of u upto which the 
structures remain pulsationally stable]. The first maxima in mass, among the equilibrium 
sequences of masses, is reached when the ratio of central pressure to central energy- density, 
<r[= (P /E )], approaches to a value about 0.6285. The binding-energy per baryon, a r [= 
(M r — M)/M r ] also approaches to its first maxima for this maximum stable value of mass. 
The corresponding maximum stable value of compactness ratio, u, is obtained as 0.3389 
(Table 2). Thus, the structure remains pulsationally stable upto a u value as large as 0.3389, 
that is, u < 0.3389 [notice that both, the upper bound on maximum value of u < 0.3406, and 
(r^o < 2.5946 are fully consistent with the corresponding absolute upper bounds]. Thus, 
the corollary 1 of theorem 2 is fully satisfied for this case. It follows, therefore, that the 
M — R relation provides the necessary and sufficient condition for dynamical stability of 
equilibrium masses, since these models are fully consistent with the compatibility criterion. 
This is evident from the comparison of column one and column seventh of Table 2 which 
indicates that for each value of (P /E ), the compactness ratio of the whole configuration, 
u, is always less than Uh, the compactness ratio of the corresponding homogeneous density 
sphere. Obviously, the upper bound on u < 0.3389 obtained here by using the compatibility 
criterion, turns out to be much stronger than the upper bounds on this parameter obtained 
by Lindblom (1984) and Haensel et al (1999). 
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5. Dynamical stability of polytropic configurations 



The dynamical stability of polytropic configurations (P = Kp Fl ) was first investigated by 
Tooper (1965) for polytropic index 3 < n < = 1 + 1/n) by using the variational 

method which states that a sufficient condition for the dynamical stability of a mass is that 
the right-hand side of the following equation 



uJ 1 



Jo 
Jo 

+ [ R e^ + ^ 2 [ 7 P/r 2 ] (r 2 e-^ 2 0' 2 dr 
Jo 

- [ R e (x+ ^[P' 2 /(P + E)]r 2 edr 
Jo 



+8tt ( R e^ x+ ^l 2 P(P + E)r 2 i 2 dr. (22) 
Jo 

vanishes for some chosen "trial function" £ which satisfies the boundary conditions 
£ = at r = 0, (23) 
and 

SP = -£P' - 1 Pe^ 2 [(r 2 e- v / 2 0'/r 2 } 

= at r = R, (24) 

where uj is the angular frequency of pulsation, R is the size of the configuration, and SP is the 
'Lagrangian displacement in pressure'. The prime denotes radial derivative, and the quantity 
7 = [(P + E)/P](dP/dE) = Ti (constant) for the polytropic configurations considered here. 
Tooper (1965) used the trial function of the form £ = &ir(l + air 2 + a 2 r 4 + a 3 r 6 + ...)e v l 2 , 
where ai,a 2 ,a 3 , ... are adjustable constants, in Eq.(22) and showed that for 3 < n < 1, 
the first maxima among the masses in the mass-central density (or, mass-radius) relation 
approaches at the same value of central pressure to central rest-density ratio (Po/po) where 
the squared frequency of pulsation, u 2 , also becomes zero. 
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In order to verify these results in view of the discussion of section 3, we choose the 
trial function £ = re^ 4 and employ a fourth-order Runge-Kutta method to solve Eq.(22). 
The results of this iteration are presented in Tables 3-4 for the polytropic index n = 1 
and 1.5 respectively. It is seen that the first maxima among masses in the mass-centre 
density relation is reached for the same value of (P /Po) where u 2 also approaches zero. 
This finding is in perfect agreement with those of the Tooper (1965) and together with 
theorem 2 verifies further that just the choice of the trial function £ = re"/ 4 is capable of 
providing the necessary and sufficient condition for dynamical stability of masses. 

6. Results and conclusions 

We have re- investigated the core-envelope model with stiffest equation of state [speed of 
sound equal to that of light] in the core and a polytropic equation with constant adiabatic 
index r\ = [dlnP/dlnp] in the envelope, based upon the criterion obtained by Negi and 
Durgapal (2001). We find that the condition of hydrostatic equilibrium is assured only when 
the minimum ratio of (Pb/E b ) at the core-envelope boundary reaches about 2.9201 x 10 _1 
[that is, when the value of the adiabatic index, F 1 at the core-envelope boundary reaches 
around 4.4246 as compared to the previous case of Ti = 72.4286]. Under this condition, 
the pressure, density, both of the metric parameters including their first derivatives, and 
the speed of sound are continuous at the core-envelope boundary and at the surface. The 
mass-radius diagram indicates that the configuration remains dynamically stable upto a u 
value as large as 0.3389. The corresponding central value of the local (r^o is obtained as 
2.5911. These values are fully consistent with those of the absolute values, u max = 0.3406, 
and (ri)o,max — 2.5946, compatible with the structure of general relativity, causality, and 
dynamical stability, obtained by using a causal configuration of homogeneous energy-density 
(Negi 2004a). Not just for the particular model considered in the present study, we have 
also found that the maximum value of u and the corresponding value of (r^o upto which 
the configurations remain pulsationally stable can not exceed the values of u ~ 0.34 and 
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(Ti)o ~ 2.5128 respectively, if the 'compatibility criterion' is followed and the envelope of 
the present model is replaced by the polytropic EOS (4/3) <dlnP/dlnp < 2 (Negi 2005). 

In addition to the study of two-density models considered in the present study (sec. 4), 
the dynamical stability of the polytropic configurations (sec. 5) explicitly shows that the 
compatibility criterion (Negi & Durgapal 2001) alone is capable of providing a necessary 
and sufficient condition for any regular configuration to be consistent with the structure of 
general relativity [however, the study of M — R relation of such sequences (corresponding 
to two-density models with v = c = 1 at the centre), consistent with the definition of actual 
mass (Negi 2004 b; Negi 2006) and the 'compatibility criterion' (Negi & Durgapal 2001) in 
the near future will finally settle down this issue]. 

The two-density structures are dynamically stable and gravitationally bound even for 
the value of compactness ratio, u > (1/3), thus giving a suitable model for studying the 
Ultra-compact Objects [UCOs] discussed in the literature (see, e. g., Negi & Durgapal 1999a, 
b; 2000; and references therein). The present type of studies may also find application to 
test various models of NSs based upon EOSs of dense nuclear matter, and the models of 
relativistic stellar objects like - star clusters. 
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Table 1: Properties of the causal core-envelope models, with a core given by the 
most stiff EOS, (dP/dE) = f, and the envelope is characterized by the polytropic EOS, 
(dlnP/dlnp) = T±, such that, all the parameters, P,E,u,X, and the speed of sound, 
(dP/dE) 1 ^ 2 , are continuous at the core-envelope boundary, r&, and the models satisfy the 
necessary (but not sufficient) condition for hydrostatic equilibrium. The maximum value 
of u[= (M/R) = 0.3574] for the structure is obtained [Fig. 1 of Negi & Durgapal (2000)], 
when the minimum value of the ratio of pressure to density at the core-envelope bound- 
ary, (Pb/Eb), reaches about 0.014. The maximum value of the binding-energy per baryon, 
a r [= (M r — M)/M r , where M r is the rest mass of the configuration] = 0.2441, also occurs 
for the maximum stable value of u. The subscript '0' and 'b' represent, the values of respec- 
tive quantities at the centre, and at the core-envelope boundary, zr stands for the surface 
redshift. The calculations are performed for an assigned value of the central energy- density, 
E — 1. Various values shown in the table are round off at the fourth decimal place. The 
slanted values represent the limiting case upto which the structure remains dynamically sta- 
ble. However, the model do not satisfy the necessary and sufficient condition for hydrostatic 
equilibrium, since for each assigned value of (P /E ), the compactness ratio, u, of the whole 
configuration always corresponds to a value larger than that of the compactness ratio, Uh, 
of the homogeneous density distribution (that is u > Uh)- Furthermore, the M — R relation 
does not provide the necessary and sufficient condition for dynamical stability of equilib- 
rium configurations, since the maximum stable value of u ~ 0.3574 exceeds the limiting 
value (u = 0.3406) for the corresponding centre value of (r^o ~ 2.4204 (which is, however, 
consistent with the corresponding absolute upper bound < 2.5946). 
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TABLES 



Uh 

0.1527 
0.1901 
0.2195 
0.2478 
0.2696 
0.2975 
0.3213 
0.3315 
0.3369 
0.3424 
0.3501 
0.3504 
0.3513 
0.3523 



Pq/Eq 
0.1110 
0.1562 
0.2012 
0.2564 
0.3100 
0.4000 
0.5070 
0.5661 
0.6010 
0.6410 
0.7040 
0.7070 
0.7151 
0.7238 



(ri)o 

10.009 
7.4020 
5.9702 
4.9002 
4.2258 
3.5000 
2.9724 
2.7665 
2.6639 
2.5601 
2.4204 
2.4144 
2.3984 
2.3816 



n 

0.1886 
0.2181 
0.2410 
0.2639 
0.2832 
0.3127 
0.3477 
0.3690 
0.3829 
0.4003 
0.4329 
0.4347 
0.4396 
0.4450 



E b 
0.9017 
0.8559 
0.8101 
0.7542 
0.6998 
0.6086 
0.5001 
0.4401 
0.4047 
0.3642 
0.3002 
0.2971 
0.2890 
0.2801 



R 
0.2012 
0.2282 
0.2489 
0.2712 
0.2887 
0.3176 
0.3522 
0.3733 
0.3870 
0.4043 
0.4376 
0.4395 
0.4443 
0.4497 



u 

0.1558 
0.1946 
0.2243 
0.2550 
0.2763 
0.3058 
0.3304 
0.3405 
0.3455 
0.3504 
0.3574 
0.3577 
0.3583 
0.3590 



a r 
0.1001 
0.1290 
0.1512 
0.1751 
0.1917 
0.2152 
0.2334 
0.2395 
0.2419 
0.2434 
0.2441 
0.2441 
0.2436 
0.2430 



0.2052 
0.2795 
0.3467 
0.4287 
0.4949 
0.0648 
0.7168 
0.7703 
0.7988 
0.8282 
0.8724 
0.8744 
0.8783 
0.8833 



0.3473 
0.4977 
0.6515 
0.8573 
1.0598 
1.4513 
2.0020 
2.3633 
2.6033 
2.9084 
3.4925 
3.5242 
3.6085 
3.7055 
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Table 2: Properties of the causal core-envelope models, as discussed in the present paper, 
with a core given by the most stiff EOS, (dP/dE) = 1, and the envelope is characterized 
by the polytropic EOS, (dlnP/dlnp) = Ti, such that, all the parameters, P,E,u,X, and 

1 /2 

the speed of sound, (dP/dE) 1 , are continuous at the core-envelope boundary, r b . The 
maximum value of u[= (M / R) = 0.3389] for the structure is obtained (Fig. 1), when the 
minimum value of the ratio of pressure to density at the core-envelope boundary, (Pb/Eb), 
reaches about 2.9201 x 10 _1 . The first maxima among the values of the binding-energy per 
baryon, a r [= (M r — M)/M r , where M r is the rest mass of the configuration] also occurs for 
the maximum stable value of u. The calculations are performed for an assigned value of the 
central energy-density, E = 1. Except (P /E ) and r b , all other values are round off at the 
fourth decimal place. The subscript '0' and 'b' represent, the values of respective quantities 
at the centre, and at the core-envelope boundary, zr stands for the surface redshift. The 
slanted values represent the limiting case upto which the structure remains dynamically 
stable. The model is fully compatible with the structure of general relativity, as it is seen 
that for each assigned value of (P /E ), the compactness ratio, u, of the whole configuration 
always corresponds to a value less than or equal to that of the compactness ratio, Uh, of 
the homogeneous density distribution (that is u < Uh). The maximum stable value of 
u ~ 0.3389(< 0.3406) and the corresponding central value of (ri) ^ 2.5911(< 2.5946) 
indicate that the M—R relation provides the necessary and sufficient condition for dynamical 
stability of equilibrium configurations. 
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Uh 

(10°) 
0.2628 
0.2652 
0.2846 
0.2976 
0.3072 
0.3112 
0.3188 
0.3261 
0.3274 
0.3306 
0.3355 
0.3394 
0.3407 
0.3409 
0.3411 
0.3438 
0.3497 
0.3553 



(10°) 
0.29202 
0.29800 
0.35452 
0.39993 
0.43923 
0.45699 
0.49397 
0.53331 
0.54099 
0.55999 
0.59146 
0.61902 
0.62850 
0.62936 
0.63148 
0.65148 
0.69998 
0.75294 



(Ti)o 
(10°) 
4.4244 
4.3557 
3.8207 
3.5005 
3.2767 
3.1882 
3.0244 
2.8751 
2.8485 
2.7857 
2.6907 
2.6154 
2.5911 
2.5889 
2.5836 
2.5350 
2.4286 
2.3281 



n 
(io- 3 ) 

2.68 
41.25 
128.53 
166.04 
191.90 
202.42 
222.83 
243.11 
246.97 
256.44 
272.04 
285.84 
290.65 
291.09 
292.18 
302.55 
329.29 
362.83 



E b 
(10°) 
0.9999 
0.9915 
0.9117 
0.8476 
0.7920 
0.7670 
0.7147 
0.6592 
0.6483 
0.6215 
0.5770 
0.5381 
0.5247 
0.5235 
0.5205 
0.4922 
0.4238 
0.3490 



R 

(10°) 
0.3117 
0.3132 
0.3275 
0.3400 
0.3509 
0.3564 
0.3678 
0.3801 
0.3830 
0.3894 
0.4020 
0.4128 
0.4176 
0.4171 
0.4182 
0.4276 
0.4533 
0.4888 



u 

(10°) 
0.2628 
0.2649 
0.2837 
0.2965 
0.3063 
0.3102 
0.3178 
0.3252 
0.3265 
0.3296 
0.3341 
0.3380 
0.3389 
0.3394 
0.3396 
0.3419 
0.3467 
0.3505 



(10°) 
0.1698 
0.1714 
0.1870 
0.1984 
0.2069 
0.2105 
0.2170 
0.2224 
0.2235 
0.2256 
0.2289 
0.2307 
0.2314 
0.2311 
0.2313 
0.2322 
0.2328 
0.2306 



zr 
(10°) 
0.4520 
0.4584 
0.5205 
0.5676 
0.6066 
0.6232 
0.6568 
0.6913 
0.6974 
0.7130 
0.7361 
0.7568 
0.7618 
0.7642 
0.7654 
0.7783 
0.8061 
0.8290 



(10°) 
0.9615 
0.9831 
1.2026 
1.3943 
1.5738 
1.6588 
1.8467 
2.0657 
2.1101 
2.2255 
2.4265 
2.6215 
2.6888 
2.6989 
2.7146 
2.8710 
3.2991 
3.8720 
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Table 3: Properties of the models characterized by the pure polytropic EOS, 
(dlnP/dlnp) = T 1 , for I\ = 5/3(n = 1.5). The stability of the models is judged by the 
variational method for the choice of the trial function £ = re u ^ in Eq.(22), as well as the 
mass-central density relation [equivalent to dimensionless mass (M/M*) vs. (P /E ) or 
(Po/po) ratio; where M* = (n + l) 3/2 (P /Po) n/2 /( 47r Po) 1/2 ]- It is apparently seen that the 
configurations become dynamically unstable beyond the maximum mass where the squared 
angular frequency of pulsation, u 2 , also approaches zero. Thus, it follows that the mass- 
central density relation provides the necessary and sufficient condition for dynamical sta- 
bility of equilibrium configurations. The maximum value of the binding-energy per baryon, 
a r [= (M r — M)/M r , where M r is the rest mass of the configuration] also occurs for the 
maximum value of mass. All values are round off at the fourth decimal place. The subscript 
'0' represents, the values of respective quantities at the centre. The slanted values represent 
the limiting case upto which the structure remains dynamically stable. The model is fully 
compatible with the structure of general relativity, as it is seen that for each assigned value 
of (Pq/Eq), the compactness ratio, u, of the whole configuration always corresponds to a 
value less than that of the compactness ratio, Uh, of the homogeneous density distribution 
(that is u < Uh). 
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0.0465 
0.0653 
0.0805 
0.1026 
0.1158 
0.1239 
0.1300 
0.1523 
0.1798 
0.2696 
0.3279 
0.3705 



Po/Po 
0.0500 
0.0724 
0.0916 
0.1213 
0.1402 
0.1522 
0.1615 
0.1974 
0.2462 
0.4526 
0.6453 
0.8340 



M/M* 
0.1907 
0.2153 
0.2272 
0.2359 
0.2378 
0.2380 
0.2377 
0.2346 
0.2273 
0.1924 
0.1684 
0.1523 



0.0277 
0.0343 
0.0381 
0.0414 
0.0422 
0.0424 
0.0422 
0.0405 
0.0357 
0.0021 
-0.0324 
-0.0624 



LO 2 /E 

1.1085 
0.8047 
0.5734 
0.2669 
0.0979 
0.0000 
-0.0715 
-0.3130 
-0.5694 
-1.0936 
-1.1804 
-1.1146 



u 

0.0698 
0.0904 
0.1046 
0.1221 
0.1308 
0.1356 
0.1391 
0.1501 
0.1604 
0.1750 
0.1716 
0.1635 



Uh 

0.0783 
0.1032 
0.1213 
0.1446 
0.1571 
0.1643 
0.1695 
0.1872 
0.2061 
0.2537 
0.2760 
0.2894 
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Table 4: Properties of the models characterized by the pure polytropic EOS, 
(dlnP/dlnp) = T 1 , for r x = 2(n = 1). The stability of the models is judged by the variational 
method for the choice of the trial function £ = re"^ 4 in Eq.(22), as well as the mass-central 
density relation [equivalent to dimensionless mass (M/M*) vs. (P /E ) or (P /po) ratio; 
where M* = (n + 1) 3 / 2 (P / Po) n ^ 2 /(^Po) 1 ^ 2 ]- It is apparently seen that the configurations 
become dynamically unstable beyond the maximum mass where the squared angular fre- 
quency of pulsation, u 2 , also approaches zero. Thus, it follows that the mass-central density 
relation provides the necessary and sufficient condition for dynamical stability of equilibrium 
configurations. The maximum value of the binding-energy per baryon, a r [= (M r — M)/M r , 
where M r is the rest mass of the configuration] also occurs for the maximum value of mass. 
All values are round off at the fourth decimal place. The subscript '0' represents, the values 
of respective quantities at the centre. The slanted values represent the limiting case upto 
which the structure remains dynamically stable. The model is fully compatible with the 
structure of general relativity, as it is seen that for each assigned value of (P /E ), the com- 
pactness ratio, u, of the whole configuration always corresponds to a value less than that of 
the compactness ratio, Uh, of the homogeneous density distribution (that is u < u h ). 
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0.1114 
0.1252 
0.1579 
0.2203 
0.2309 
0.2460 
0.2548 
0.2800 
0.3266 
0.3948 
0.4430 
0.4917 



Po/po 
0.1254 
0.1431 
0.1875 
0.2826 
0.3002 
0.3263 
0.3420 
0.3889 
0.4850 
0.6523 
0.7952 
0.9674 



M/M* 
0.1752 
0.1830 
0.1958 
0.2055 
0.2059 
0.2060 
0.2058 
0.2045 
0.2000 
0.1902 
0.1820 
0.1733 



0.0699 
0.0744 
0.0824 
0.0895 
0.0898 
0.0899 
0.0897 
0.0885 
0.0835 
0.0704 
0.0576 
0.0418 



U 2 /Eq 

2.2225 
1.9497 
1.3440 
0.3512 
0.2030 
0.0000 
-0.1126 
-0.4151 
-0.8980 
-1.4356 
-1.7007 
-1.8762 



u 

0.1452 
0.1559 
0.1773 
0.2077 
0.2115 
0.2167 
0.2194 
0.2265 
0.2365 
0.2458 
0.2486 
0.2490 



Uh 

0.1530 
0.1655 
0.1913 
0.2301 
0.2356 
0.2430 
0.2471 
0.2580 
0.2755 
0.2961 
0.3080 
0.3184 
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FIGURES 

Fig. 1. Mass- Radius diagram of the model corresponding to Table 2, for an assigned value of 
E = Ef, = 2 x 10 14 g cm~ 3 at the core-envelope boundary r^, such that the compactness ratio, u, 
of the whole configuration always turns out less than or equal to the compactness ratio , u^, of 
homogeneous density sphere. This requirement is fulfilled only when the ratio of pressure to density 
(Pb/Eb) at rb reaches to a minimum value about 2.9201 x 1CT 1 . The pressure, energy-density, u, A, 
and the speed of sound, (dP/dE) 1 ^ 2 are continuous at the core-envelope boundary. 
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